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Abstract. Topological phases can not only be protected by internal symmetries (e.g., 
time-reversal symmetry), but also by crystalline symmetries, such as reflection or rotation 
symmetry. Recently a complete topological classification of reflection symmetric insulators, 
superconductors, nodal semimetals, and nodal superconductors has been established. In this 
article, after a brief review of the classification of reflection-symmetry-protected semimetals 
and nodal superconductors, we discuss an example of a three-dimensional topological Dirac 
semimetal, which exhibits time-reversal symmetry as well as reflection and rotation symmetries. 
We compute the surface state spectrum of this Dirac semimetal and identify the crystal lattice 
symmetries that lead to the protection of the surface states. We discuss the implications of 
our findings for the stability of the Fermi arc surface states of the Dirac material NaaBi. Our 
analysis suggests that the Fermi arc of NasBi is gapped except at time-reversal invariant surface 
momenta, which is in agreement with recent photoemission measurements. 


1. Introduction 

Finding new topological materials beyond the Bi-based topological insulators has become 
an important goal in solid state research. In the past few years, it has been realized 
that not only internal symmetries P 13 13 H E], such as time-reversal, but also crystal 
symmetries, for example crystal reflection or rotation, can lead to the protection of topological 
states is m 0 ei mi u m o nu ns uni n?]- In fact, a complete classification of reflection- 
symmetry-protected insulators and fully gapped superconductors was recently presented in 
Refs. 03 U) IU3- The topological character of these crystalline topological insulators and 
superconductors gives rise to gapless states at the boundary. Importantly, since the surface of 
a material has lower crystal symmetry than its bulk, not all surfaces of a crystalline topological 
insulator (or superconductor) exhibit gapless boundary states. Only those surfaces that are 
invariant under the crystal symmetry can support topological surface states. Experimentally, 
SnTe has been identified to be a crystalline topological insulator with Dirac cone surface 


states, that are protected against gap opening by reflection symmetry together with time- 
reversal symmetry mmmm- Other candidate materials for reflection-symmetry-protected 
topological insulators are the anti-perovskites CaaPbO and SraPbO [2211231122] • 

The concept of topological band theory can also be applied to nodal systems, such as 
semimetals and nodal superconductors. The stability of the Fermi surfaces (superconducting 
nodes) of these materials is ensured by the conservation of a topological invariant which is 
defined in terms of an integral along a contour enclosing (encircling) the Fermi point (Fermi line). 
Topological nodal systems can be protected by global symmetries [251 ESI 271 281 [291 [3Q1 ED [32] 
(i.e., internal symmetries, such as time-reversal) as well as crystal lattice symmetries, or a 
combination of the two [33l E2 [351136]. The nontrivial topology of the electronic band structure 
of these gapless materials manifests itself at the boundary in terms of protected surface states. 
Depending on the symmetry properties and the dimension of the Fermi surface, these surface 
states form Dirac cones, flat bands, or Fermi arc states |29l . Recently, several materials have 
been experimentally identified to be crystal-symmetry-protected topological semimetals. Among 
them are the Dirac materials CdsAs 2 [33 |38l [391 Hol SD S21H3] and NasBi mmmm, whose 
Dirac points are protected by rotation symmetry. The Fermi arc states of NasBi have recently 
been observed in angle-resolved photoemission measurements [26J |47]. 

In this paper, we first present in Sec. [2] a short review of the topological classification of 
reflection symmetric semimetals and nodal superconductors m- For brevity, we focus on the 
case where the Fermi surfaces (superconducting nodes) are located within the reflection plane. 
Second, we discuss in Sec. [3] a four-banc! tight-binding model on the cubic lattice describing a 
three-dimensional Dirac semimetal, which is protected by C 4 rotation symmetry. We discuss 
the surface states of this topological semimetal and determine the symmetries that protect the 
gapless surface states. Our tight-binding model has similar symmetry properties and topological 
features as the Dirac semimetal NasBi. Hence, from our topological analysis we can obtain 
information about the stability of the Fermi arc surface states of NasBi, see Sec. |3.4| 

2. Topological classification of reflection symmetric semimetals and nodal 
superconductors 

Before discussing a concrete example of a crystalline topological semimetal, let us briefly review 
the topological classification of reflection symmetric semimetals and nodal superconductors m- 
The topological properties of reflection symmetric semimetals (nodal superconductors) can 
be classified in a similar manner as those of reflection symmetric insulators (fully gapped 
superconductors) [151 16] [17]. The classification depends on the internal symmetry properties of 
the nodal system as well as on whether the reflection symmetry commutes or anticommutes with 
the internal symmetries. There are three fundamental internal symmetries, which act locally in 
position space, namely time-reversal symmetry (TRS), particle-hole symmetry (PHS), and chiral 
or sublattice symmetry (SLS). In momentum space, time-reversal and particle-hole symmetry 
act on the Bloch Hamiltonian H( k) as 


T -1 f7(—k)T = +if(k) 

(la) 

C~ 1 H(-'k)C = -H{ k), 

(lb) 


respectively, where T = JCUt and C = JCUc are antiunitary operators and 1C is the complex 
conjugation operator. Chiral symmetry is implemented in terms of a unitary matrix S which 
anti-commutes with the Hamiltonian, i.e., SH(k) + H(k)S = 0. 

Crystal reflection symmetry acts non-locally in position space. In the following, we consider 
a dez-dimensional Hamiltonian H( k) which is symmetric under reflection in the first direction. 


Hence, H( k) satisfies 


R~ 1 H(-k 1 ,k)R = H(ki,k), (2) 

where k = (& 2 , & 3 , • • •, kd BZ ) and i? is a unitary reflection operator. Since a phase factor can be 
absorbed in the definition of the electron creation and annihilation operators, we can assume 
that R is Hermitian, i.e., R) = R. As a consequence, the eigenvalues of R are either +1 or — 1 
and all commutation and anti-commutation relations between R and the internal symmetries 
are defined unambiguously, i.e., 

TRT- 1 = 7) T R, CRC- 1 = rj c R and SRS _1 = r] S R, (3) 

where the three indices t]t, tjc, and 17 s specify whether R commutes (+ 1 ) or anti-commutes (— 1 ) 
with the corresponding internal symmetry operator. For the symmetry classes AI, All, AIII, C, 
and D these different possibilities are labeled in Table [l] as R Vt , R Vc , and R Vs , respectively. For 
the symmetry classes that contain two internal symmetries (i.e., classes BDI, Cl, CII, and Dill), 
the different (anti-)commutation relations are denoted by R VrVc . Hence, reflection symmetries 
together with the three internal symmetries define a total of 27 symmetry classes, see Table [lj 
The classification of reflection-symmetry-protected semimetals and nodal superconductors 
depends on the symmetry class and on the codimension of the Fermi surface (superconducting 
nodes) p = dpz ~ dps, where dpz denotes the spatial dimension (i.e., the dimension of the 
Brillouin zone) and dps is the dimension of the Fermi surface. Furthermore, the classification 
depends on how the Fermi surfaces (superconducting nodes) transform under the reflection 
symmetry and the internal symmetries. In general, one can distinguish three different situations: 
(i) Each Fermi surface is left invariant by both reflection and internal symmetries, (ii) Fermi 
surfaces are invariant under reflection symmetry, but are pairwise related to each other by the 
internal symmetries, and (iii) different Fermi surfaces are pairwise related to each other by both 
reflection and in internal symmetry operations. In case (i) and case (ii) the Fermi surfaces are 
located within the refection plane, whereas in case (iii) the Fermi surfaces are positioned away 
from the reflection plane. For brevity we focus here on only case (i) and case (ii). Case (iii) has 
been discussed extensively in Refs. [3.3j and [35]. 

2.1. Fermi surfaces at high-symmetry points within mirror plane 

First we consider case (i), i.e., Fermi surfaces that are invariant under both reflection and 
internal symmetries. These Fermi surfaces are located within the reflection plane and at high- 
symmetry points in the Brillouin zone, that is, at time-reversal invariant momenta. For this case 
the classification of reflection symmetric Fermi points can be inferred from the classification of 
reflection-symmetry-protected insulators by means of a dimensional reduction procedure [ 33] . 
Namely, the surface states of reflection symmetric dez-dimensional topological insulators can be 
viewed as reflection-symmetry-protected Fermi points in dp z — 1 dimensions. From this it follows 
that the topological classification of reflection symmetric Fermi points at time-reversal invariant 
momenta with dps = 0 is obtained from the classification of reflection symmetric insulators by 
the dimensional shift dp z —> dpz ~ 1; see Table [lj This logic also works for Femi surfaces with 
dps > 0, if their stability is guaranteed by an MZ or 2MZ topological number. However, Z 2 
and AfZ 2 topological numbers ensure only the stability of Femi points, i.e., Fermi surfaces with 
dpg = 0. Derivations based on Clifford analysis and Iv theory HZ1I33] corroborate these findings. 
The classification of reflection-symmetry-protected Fermi surfaces at high symmetry points of 
the Brillouin zone is presented in Table [lj where the first row indicates the codimension p of the 
Fermi surface. We note that the classification is eight-fold periodic in p. 


Table 1 . Topological classification of reflection-symmetry-protected semimetals and nodal 
superconductors. The first and second rows specify the codimension p = d-^z — g^fs of the 
reflection symmetric Fermi surfaces (superconducting nodes) at high-symmetry points and 
away from high-symmetry points of the Brillouin zone, respectively. The third row shows the 
classification of reflection symmetric insulators and fully gapped superconductors [El ESI ■ The 
first column indicates whether the reflection operator R commutes (“R+”, or anti¬ 
commutes “/?,_”) with the global symmetries. “R^ _” and “R _|_” denote the case 

where R commutes with one of the global symmetries but anti-commutes with the other one. 
The second column lists the global symmetry classes (using the Cartan nomenclature |3|), which 
are distinguished by the presence or absence of time-reversal symmetry, particle-hole symmetry, 
and chiral symmetry. 
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Z 2 and MZ 2 invariants only protect Fermi surfaces of dimension zero (i.e., Fermi points with dps — 0) at high-symmetry points of the Brillouin zone. 
Fermi surfaces located within the mirror plane but away from high symmetry points cannot be protected by a Z 2 or MZ 2 topological number. However, 
the nodal system can exhibit zero-energy surface states that are protected by a Z 2 or MZ 2 topological number. 

For topological semimetals (nodal topological superconductors) the presence of translation symmetry is always assumed. Therefore, one does not need 
to distinguish between TZ 2 and Z 2 invariants for these gapless topological materials. 


2.2. Fermi surfaces within mirror plane but off high-symmetry points 

Second, we review case (ii), i.e., the classification of reflection symmetric Fermi surfaces that are 
located away from high-symmetry points of the Brillouin zone. These Fermi surfaces transform 
pairwise into each other by the global (i.e., internal) symmetries, which relate k —> — k. Using 
an analysis based on the minimal-Dirac-Hamiltonian method, it was show in Ref. |33] that 
only MZ and 2MZ topological numbers can ensure the stability of reflection symmetric Fermi 
surfaces off high-symmetry points. Z 2 and MZ 2 invariants, on the other hand, do not give rise to 
stable Fermi surfaces. However, as illustrated in Sec. [3j Z 2 or M Z 2 invariants may nevertheless 
lead to protected zero-energy surface states at time-reversal-invariant momenta of the surface 
Brillouin zone. The classification of reflection symmetric Fermi surfaces located away from high- 
symmetry points is summarized in Table [l] As above, we note that due to Bott’s periodicity 
theorem |48j , the classification scheme exhibits an eight-fold periodicity in p. We observe that 
the classification of reflection-symmetric Fermi surfaces located away from high symmetry points 












with codimension p is related to the classification of reflection-symmetric insulators with spatial 
dimension dez = p — 1 • 


3. Topological Dirac semimetal with reflection and rotation symmetry 

In this section we discuss a tight-binding model describing a topological semimetal which exhibits 
time-reversal symmetry as well as reflection and rotation symmetries. This model is inspired by 
the Dirac material NasBi, whose three-dimensional bulk Dirac cones are protected by reflection 
symmetry. Using angle-resolved photoemission experiments, topological surface modes have 
recently been observed on the (100) surface of NasBi [ 3T]| . 


3.1. Model Hamiltonian 

We consider the following cubic-lattice Hamiltonian describing a four-band semimetal with two 
Dirac points 



[^?+z ( r2 so+iT x s x )Cg-\-Ct + ~ 


( T z s 0 +iT x Sy)C £ +Ct + .T z soC £ +mCtT z soC £ +h.c., (4a) 


with the spinor Cg = [ci-|-(af), c\±(x), C 2 f(x), C 2 i(x)], where c nc7 (x) denotes the electron 

annihilation operator with orbital index n and spin a. The two sets of Pauli matrices r Q and 
s a operate in spin space and orbital space, respectively. In momentum space Hamiltonian (??) 
takes diagonal form % = Sk Ck-ff(k)C k , with C k the Fourier-transformed spinor, and 

H (k) = sin k x T x s x + sin k y T x s y + (cos k x + cos k y + cos k z — m)r 2 so- (4b) 


For concreteness we set m = 2. With this choice, the bulk Dirac cones of H( k) are positioned at 
k± = (0,0, ±vr/2), i.e., they are located on the k z axis just as in the Dirac material NasBi. The 
cubic-lattice Hamiltonian (4b) possesses a rotation symmetry of order four about the z axis. 


This C\ rotation symmetry acts on H( k) as 


R C4 H( ky^k x ^k z )R(j^ — H(k x , k y , k z ), 


(5) 


where Rc 4 = t 2 (so + is z )/\/ 2 denotes the C 4 rotation operator. The Dirac semimetal © also 
preserves time reversal symmetry, Eq. (??), with the time-reversal operator T = tqs v 1C. Since 
T 2 = —1, Hamiltonian Q belongs to symmetry class AIL The codimension of the Dirac cones 
of H( k), which are located away from the high symmetry points in the Brillouin zone, is p = 3. 
Hence, according to the classification of Ref. m (see Table 1 in Ref. E33), the Dirac points are 
not protected by time-reversal symmetry, even though a binary Z 2 invariant can be defined for 
this symmetry class. (As we will show below, this Z 2 invariant leads to gapless surface states 
at time-reversal-invariant momenta of the surface Brillouin zone.) Indeed, we find that the 
bulk Dirac cones can be gapped out by the mass term smk z T y so, which preserves time reversal 
symmetry. However, the mass term sin k z T y So breaks the C 4 rotation symmetry ([ 5 ]). In general, 
we find that the only possible gap opening terms are fi(k z )T y so and f r 2 (k z )T x s z , since these 
are the terms that anti-commute with H( k). Here, fi(k z ) and f 2 (k z ) represent k z dependent 
masses. We find that these two gap opening terms break the C4 rotation symmetry ([ 5 ]), since 
they anti-commute with Rc 4 - As a consequence, the two Dirac cones of Hamiltonian Q are 
protected against gap opening by the C4 rotation symmetry ([ 5 ]). 


3.2. Surface states 

As discussed above and in Ref. [33], class All Dirac points with codimension p = 3 that are 
located away from high-symmetry points of the Brillouin zone are unstable in the absence of 





Figure 1. (a) Surface band structure of the Dirac semimetal Q for a (100) slab as a function 

of surface momentum k y with fixed k z = 7t/ 2. (b) Same as panel (a) but for k z = 0. The 
gapless surface state at k y = 0 is protected by a Z 2 invariant, (c) Band structure at the (100) 
surface of the Dirac semimetal Q as a function of k y with k z = 7r/4. In the absence of reflection 
symmetry ([o]) and chiral symmetry (JT]) the surface state (red trace) is gapped, (d) Surface band 
structure of the Dirac semimetal Q (i.e., with reflection symmetry and chiral symmetry) as a 
function of k y with k z = n/ 4. 


crystal lattice symmetries, even though a binary Z 2 index can be defined for this symmetry 
class. This Z 2 invariant does not guarantee the stability of Dirac points, however it gives rise 
to protected gapless surface states at time-reversal invariant momenta of the surface Brillouin 
zone. For the Dirac semimetal (0 these gapless modes appear at ( k y ,k z ) = (0,0) of the (100) 
surface Brillouin zone, see Fig.lHb). For other surface momenta, the surface states are in 
general gapped. Similarly, rotation symmetry 0 can protect surface states only at time-reversal- 
invariant momenta of the surface Brillouin zone, but not away from these points. 

To illustrate this, we consider Hamiltonian © with periodic boundary conditions along the 
z-axis and open boundary conditions along the x- and y-directions. We observe that model 
Hamiltonian @ exhibits besides rotation symmetry 0 an accidental reflection symmetry 

R y H{k x , ky, k z )R y — H(k x , k y , k z ), (6) 

with R y = T z s y , and an accidental chiral symmetry 

H(k)S = -SH( k), (7) 


with S = t x s z . In order to break these accidental symmetries we need to add to the Hamiltonian 
the term +g sink z T x s z on the (100) and (100) surfaces, and the term — g sin k z T x s z on the (010) 
and (010) surfaces. That is, we consider 


H(k z ) = H(k z )+g 
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Cl ± sin k z r x s z 
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X±,x\_ 
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( 8 ) 


where H{k~) is obtained from Eq. (??) by Fourier transforming the z-coordinate and the sum 
in Eq. ([8]) is over a collection of one-dimensional chains which are oriented along the z-axis 
and labeled by x±. Hamiltonian ([8]) is symmetric under time-reversal symmetry and rotation 
symmetry ([ 5 ]), but breaks the accidental symmetries ([6]) and Fig.g c) shows the (100) 

surface spectrum of Hamiltonian ([8]) with g = 0.2 as a function of surface momentum k y for 
k z = 7t/ 4. Fig. [2]ja) displays the (100) surface spectrum as a function of both surface momenta 
k y and k z . We observe that the surface state is gapped, except at ( k y ,k z ) = (0,0), where it is 
protected by time-reversal symmetry (i.e., the aforementioned Z 2 invariant). 
























Figure 2. (a) Surface spectrum of the Dirac semimetal ([8]) for the (100) face as a function of 

surface momenta k y and k z . The surface states and bulk states are colored in green and gray, 
respectively. In the absence of reflection symmetry ([6]) and chiral symmetry 0 the surface state 
is gapped except at (k y ,k z ) = (0,0). (b) Surface spectrum of the Dirac semimetal i.e., in 

the presence of reflection symmetry and chiral symmetry. The gaplessness of the Fermi arc is 
guaranteed by the mirror winding number Q. 


On the other hand, if we consider a semimetal which possesses also symmetries ([h]) and ([t]) 
[i.e., Hamiltonian Q instead of Hamiltonian ([8])], then the surface states are gapless along a 
one-dimensional arc that connects the two projected Dirac nodes k± = (0, 0, ±7 t/2 ), see Fig.[l](d) 
and Fig. [2j[b) . Since the reflection operator R y anti-commutes with the time-reversal operator 
T = Tos y 1C and with the effective particle-hole operator C = T~ l S = ir x s x IC , Hamiltonian Q 
can be viewed as a member of symmetry class Dill with R-- in Table [lj Hence, according to the 
classification scheme of Sec. 2.2 the Dirac nodes of Eq. (0 are protected by an M Z invariant, 
i.e., a mirror winding number. Moreover, Hamiltonian (H exhibits Fermi arc surface states, 
which are protected by this mirror winding number. This mirror invariant is defined in terms 
of a one-dimensional integral along a contour that lies within the mirror plane k y = 0. Within 
the k y = 0 mirror plane the Hamiltonian can be block diagonalized with respect to R y . For the 
block with mirror eigenvalue R y = +1 the mirror winding number takes the form 


Ur = — 
c 2vr 


Tr (q*dq), with q = 


sin k x — iM(k x ) 
sin 2 k x + M 2 {k x ) 


(9) 


and M(k x ) = cos k x + cos k z — 1, where C denotes a contour that lies within the mirror plane. 
Choosing the contour C to be parallel to the k x -axis with k y = 0 and k z a fixed parameter we 
find that 


n + (k z ) 


-1, \k z \ < 7t/2 

0, \k z \ > 7t/2 


( 10 ) 


This indicates that there exists a gapless Fermi arc state at k y = 0 on the (100) surface, which 
is in agreement with the spectrum shown in Fig. [2](b). 


3.3. The role of translation symmetry 

While each Dirac point of a topological semimetal is stable against deformations that are local 
in momentum space, the Dirac points are not protected against commensurate perturbations, 
such as charge-density wave modulations, which connect Dirac points with opposite momenta. 




However, these deformations are forbidden, if we impose besides internal and crystal point-group 
symmetries, also translation symmetry. 

To illustrate the instability of the Dirac points against charge-density wave perturbations, let 
us consider a low-energy description of semimetal Q. Expanding Eq. Q around the two Dirac 
points k-t = (0, 0, ± 7 r/ 2 ), we obtain the following low-energy effective Hamiltonian 


H eS 


E ( C K+ cL)( h Z 

kx fey j&A 



( 11 ) 


where Kj- = ( k x , k y , k/±±n/2) and h± = k x t x s x + k y t x s y pk^T z so• Hamiltonian H c g- is symmetric 
under the C 4 rotation symmetry which transforms (7 k ± —> t z (sq + is z )C^ /\/ 2 . and the 
time reversal symmetry (??), which transforms Ck ± —>• ToSy/C(7_K ± • We find that the two Dirac 
cones can be gapped out by the charge-density wave modulation 


f= E K 

kx iky ,/ca 




T z So\ f Ck+A 

0 J [c K _J 


( 12 ) 


which preserves rotation symmetry but breaks translation symmetry. We conclude that 
in the absence of translation symmetry, the Dirac cones of semimetals are unstable against 
charge-density wave-type perturbations. 


3-4- Implications for Na^Bi 

Albeit NasBi has hexagonal Ph^/mmc crystal structure rather than cubic one, its topological and 
symmetry properties are similar to the model Hamiltonian Q . Instead of C 4 rotation symmetry, 
NaaBi exhibits a C 3 rotation symmetry, which protects the two gapless bulk Dirac cones. Using 
a downfolding procedure, it is possible to derive from DFT band structure calculations a low- 
energy hexagonal-lattice tight-binding model for NaaBi. This tight-binding model is expected 
to exhibit similar symmetry properties and the same topological features as our tight-binding 
Hamiltonian ([ 8 ]). Hence, this suggests that NasBi supports Fermi arc surface states, which are 
gapped except at time-reversal invariant points of the surface Brillouin zone. This expectation 
is confirmed by recent angle-resolved photoemission experiments |%S1 ITT] . 

4. Summary and discussion 

In this paper, we have reviewed the topological classification of reflection-symmetry-protected 
semimetals and nodal superconductors. This classification scheme depends on (i) the internal 
symmetry properties of the Hamiltonian, (ii) whether the reflection symmetry commutes 
or anti-commutes with the internal symmetries, (iii) the codimension of the Fermi surface 
(superconducting node) p = dsz — ^FS; and (iv) the transformation properties of the Fermi 
surface (superconducting node) under the action of the reflection and internal symmetries. 
Reflection symmetry together with the internal (i.e., non-spatial) symmetries define a total 
of 27 symmetry classes. The result of the classification scheme is presented in Table [l] The 
stability of the Fermi surface is ensured by the conservation of a topological invariant that is 
defined as an integral along a contour surrounding the Fermi surface. These topological numbers 
can be Chern or winding numbers (indicated by “Z” in Table [I]), binary invariants (indicated 
by “Z 2 ”), mirror Chern or mirror winding numbers (indicated by “M Z”), or mirror binary 
invariants (indicated by “M Z 2 ”). 

To illustrate the usefulness of the classification Table [lj we have discussed in Sec. [3] an 
example of a topological semimetal with time-reversal symmetry and rotation symmetry. This 
model Hamiltonian represents a cubic lattice version of the hexagonal Dirac material NasBi. We 


determined the surface state spectrum of this Dirac model system and discussed its topological 
properties. We found that the (100) surface supports Fermi arc surface states, which are gapped 
except at the T point of the surface Brillouin zone, see Fig. [2](a). This is in agreement with 
recent angular-resolved photoemission experiments on the (100) surface of NasBi 1(11 T7]. 
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